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We present a calculation of the modulation in the Local Density Of electronic States (LDOS) 
caused by an impurity in graphene in the presence of external magnetic field. We focus on the 
spatial Fourier Transform (FT) of this modulation around the impurity The FT due to the low 
energy quasiparticles are found to be nonzero over the reciprocal lattice corresponding to graphene. 
At these lattice spots the FT exhibits well-defined features at wavevectors that are multiples of the 
inverse cyclotron orbit diameter (see Figure |2| and is cut off at the wavevector corresponding to the 
energy of observation. Scanning Tunneling Spectroscopy (STS) on graphene and the energy-resolved 
FT fingerprint obtained therefrom may be used to observe the quasiparticle interference of Dirac 
particles in graphene in the presence of magnetic field. 
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Graphene is a monatomic layer of carbon atoms ar- 
ranged in a hexagonal lattice. It was first isolated by the 
mechanical exfoliation of graphite in 2004 [1 . The elec- 
tronic band structure of graphene is characterized by two 
points K and K* (at wavevectors ±KT) [2 in the recipro- 
cal space where the valence and conduction bands touch 
each other. The gapless low energy excitations that exist 
at those points can be described by theories of massless 
Dirac quasiparticles with opposite chirality in (2+1) di- 
mensions [3 . An interesting consequence of such a band 
structure is the formation of Landau levels (in a perpen- 
dicular magnetic field) whose energies vary as the square 
root of the Landau level index as well as that of the mag- 
nitude of the perpendicular magnetic field [4j |5] . The un- 
conventional Quantum Hall Effect seen in transport mea- 
surements in graphene is another profound physical con- 
sequence of the Dirac nature of these quasiparticles [6j [7]. 
Till date, the only evidence of the Landau quantization 
of Dirac particles has come from bulk transport measure- 
ments. Alternately, one might look at the evidence from 
Landau level spectroscopy using the Scanning Tunneling 
Microscope (STM). In addition to the real space imaging 
of Landau levels local spectroscopic tools can be a very 
sensitive probe of QuasiParticle Interference (QPI) that 
often reveals details about the underlying band structure 
and the quasiparticle wavefunctions[HJ|9]. Applications of 
these ideas to graphene are natural and promising. Ex- 
periments are currently in progress that probe the sig- 
natures of QPI in graphene in the presence of a mag- 
netic field. In this paper we focus on QPI in magnetic 
field in graphene. We use the theory of non-interacting 
Dirac quasiparticles in a magnetic field and calculate the 
change in the electronic LDOS in response to weak im- 
purities. We find that i) the LDOS FT displays charac- 
teristic rings whose size is set by the inverse cyclotron 
diameter d cyc - indeed, in the limit of a strong magnetic 
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FIG. 1: Comparing the Green's functions when B — 
(dashed gray) and B / (continuous black), at an energy 
E ohs corresponding to the Landau level index ull = 8. The 
distance propagated is measured in units of the classical cy- 
clotron orbit diameter d cyc = 2\E\/(evB) = 2\/2ull£b. At 
small distances these oscillate together at the wavevector 
E ohs /(hv) but after ull/2 oscillations the green's function 
for B y£ decays exponentially since the particle 'turns' in its 
cyclotron orbit and cannot propagate further than d cyc - 

field when effects of disorder and line broadening are sec- 
ondary the main feature of quasiparticle motion will be 
the cyclotron orbits - and ii) these rings will form a lat- 
tice in Fourier space that is the same as the graphene 
reciprocal lattice. Also, depending on the detailed impu- 
rity potential structure this FT could have additional an- 
gular dependence in /c-space determined by off-diagonal 
(sublattice mixing terms) in impurity scattering matrix. 
These ring-like signatures could be observed in STM ex- 
periments. 

In this work we have followed the lattice-related con- 
ventions used in [2|. The hamiltonian near the if*-point 
is related to that near the if-point by a parity trans- 
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FIG. 2: Comparing the angular averages of the spatial FT and power spectrum of LDOS modulations around a short-ranged 
impurity potential Y(k) oc I, for the cases when the magnetic field is zero and when it's nonzero (and the nearest Landau 
level has an index till — 8). The center figure (C) shows part of the reciprocal lattice formed by regions in /c-space where the 
Fourier transform may be nonzero. The green and red 'spots' arise from K — > K' scattering and vice versa respectively. The 
grey spots arise from intravalley scattering. One green region is enlarged to show the angle- averaged Power Spectrum on a scale 
where the oscillations are better resolved (the density maps have edges of length 6E ohs /(hv)), for the B ^ (Rl) and B = 
(LI) cases. Below these are the corresponding variation of the Fourier transforms with k — the deviation from the if -point, 
for the B / (R2 - see ^ for parameters) and B = (L2) cases. All four plots were made using MathematicafTO] . 



formation of the lattice: H(—K + k) = a x H(K — k)a x . 
Since the low energy theory describing free excitations 
near the K/K* points obey a (2+l)-dimensional Dirac 
theory [3l 0], we can again use the parity operator <j z for 
the dirac fields within a given valley to relate the station- 
ary eigenstates of and the contributions to the propaga- 
tor/Green's function from the two valleys (modulo the 
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factors) : 



G-K = 0-yGjKO-\ = -CTyGKO-y 



(1) 



The total low-energy electronic green's function will fi- 
nally be given by: 
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(2) 



When a finite perpendicular magnetic field Bz is present, 
in a convenient gauge choice the energy eigenstates at the 
if-point are given by (for B > 0): 
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where a n = sgn(n)(l - S n0 ), £ = yjK/(eB), j n = 2- <5 n>0 , 
co c = y/2v/£, E n ^ = o- n yf\a\hjj c and n G Z. The 



(/> n 's are the orthonormal eigenfunctions of the ID Sim- 
ple Harmonic Oscillator [T3] (take (/>_i = 0). The fc's 
are consistent with periodic boundary conditions in the 
x-direction. Since the SHO eigenfunctions need to be con- 
fined inside the sample, we end up with a degeneracy of 
M = L x L y /(27rl 2 ) per Landau level. 

The green's function may be calculated as: 
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where p and are the modulus and the argument respec- 
tively of the complex number g = (x' — x) + i(y' — y) 
and q = (x r — y 1 —y)] the argument of the (associated) 
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Laguerre polynomials, denoted by L, is ^ and we have 
defined the matrices M and N for later reference. 

We note here that we haven't found this convenient 
elementary result in the literature till date. 

The K* point eigenfunctions as well as the Green's 
functions are related to those at the K point by the afore- 
said relations 0. 

It is instructive to compare the zero field and finite field 
free electronic green's functions. A comparative plot of 
the angular average of one of their components is shown 
in Figure [I] plotted against the spatial separation as a 
fraction of the cyclotron diameter d cyc = 2\E\/(evB) = 
2^2n~ZLl B . 

We now consider the case when an impurity potential 
Y(r) (whose spatial fourier transform is given by Y(q)) 
is present in graphene. We consider the general form of 
Y(r) in what follows — i.e, V(r) is a general hermitian 
2x2 matrix function of r (Y(r) = Y^(r)). 

Scanning Tunneling Microscopes (STMs) give a signal 
corresponding to the local value of the spectral function 
[TT] . We can use the Green's function derived above to 
obtain the spatial FT of the change in the LDOS (given 
by the change in the spectral function) to the linear order 
in the impurity potential strength as follows: 
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SA(k, u) = -23Tr X(fc, u + zr/)V(fe) 



(5) 



In the above, X(fc, z) is the fourier transform (w.r.t q) of 
X(g,z)=N(-g,z)N(g,z) 

= N K (-g, z)N K (g, z) + N~ k (-q, z)N~ K (g, z) 

+ N K (-Q, Z)N~ K (Q, Z) + n~ K (-Q, z)n K (Q, Z) (6) 

From Q and Q we see that N ±K (g, z) possess the pref- 
actors e ±lK ' e . We thus deduce that scattering by the 



impurity potential V will yield spatial LDOS oscillations 
around the wavevectors (terms in the first row of (J6|) 
and ±2K (second row of (|6|) as well as those joined to 
these by reciprocal lattice vectors, due to intra and in- 
tervalley scattering respectively. The resulting lattice is 
identical to the reciprocal lattice (see Figure |2^C)). 

Since scattering around the zero wavevector can also 
arise from many slowly varying unknown environmental 
potentials, we expect that LDOS oscillations around the 
wavevectors ±2K near isolated atomically sharp defects 
will better reproduce the LDOS profiles that our theory 
predicts and for this reason we'll focus on explaining how 
to calculate the features around ±2K. To do this we need 
to isolate in ([5| the part due to intervalley scattering, 
which amounts to using the terms in the last line of (|6| 
that we shall refer to as e T<2lK ' e X T respectively. 

The Landau levels have been assumed to be sharp in 
the treatment so far and so direct evaluation of Q and 
the subsequent calculations will yield a sum of delta func- 
tions in energy. To be able to resolve the spatial func- 
tional forms and to reflect realistic experimental condi- 
tions we can either assume that the Landau levels are 
broadened or that the STM has a finite detection win- 
dow. We have chosen to take a gaussian detection win- 
dow with width ri!4l: 
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It is now possible to write down the LDOS as a series ex- 
pansion in ^£ , where AE is of the order of the difference 
between the energy levels that are incorporated into the 
calculation. To see this, we note that upon substituting 
the expression (|4| of the green's function we come across 
sums of the following structure (g represents the gaussian 
in ([7]); f mn is proportional to Fourier transforms of the 
form Jd 2 qTr[M^ K (q-k)M^ K (q)Y(T^K + k)} that sat- 
isfy the condition ^/ mn = —^sfnm when V is invariant 
under spatial inversion): 
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The 'diagonal' term involving f nn above (that corre- 
sponds to the particle ejected by the STM tip remaining 
in the same Landau level on both legs of its journey before 
and after scattering off the defect) gives the main contri- 



bution and the other members in the sum are supressed 
by the aforesaid factors of The numerical calcula- 
tions that we subsequently perform are taking only the 
first few terms of this series into account and work well 



4 



for large magnetic fields when AE oc \[B ^> T. 

From ([5| and subsequent discussions we find that 
we can write 5A ohs in the form 5A(±2K + k,cu) ohs = 
TrpD>±(k, u)V(±2K + fc)], where B± = i(X± - X^) . We 
can make the following general comments regarding the 
functional dependance of the components of B±(fc, uu) as 
a function of k. Let ull(^) denote the Landau level 
index corresponding to the Landau level nearest to the 
energy of observation E ohs = hw — it is thus the in- 
teger closest to sgnuo(ou /ou c ) 2 ). The 'diagonal' term in 
(|8| that is the most important contribution then corre- 
sponds to n = till(^)- From the definition of D and 
using Q we see that D consists of products of two oscil- 
latory functions (like those shown in Figure [T]) . We thus 
expect spatial oscillation scales set by the wavevectors 
2ir/ (2^/2 \ull\£) and twice of \oo\/v to appear in D(r). 
Our calculation confirms this expectation — we find that 
D(fc) displays a set of about \ull\ oscillatory peaks start- 
ing at k = and separated by a period Ak (see below); 
it then decays rapidly after a maximum wavevector /c max , 
where 

t> y \nLL\ d cyc v i 

The off-diagonal elements in Q possess an angular 
dependence and for this reason H)(±2K + fc) exhibits si- 
nusoidal oscillations in Ok and 20 ^ for a given fc, Ok being 
the orientation angle of fc with respect to the direction 
of K. We find that when intravalley scattering is consid- 
ered, only the off-diagonal components of V give rise to 
^-oscillations while in the case of intervalley scattering, 
the diagonal components of V can, in addition, lead to 
2#/ c -oscillations. 

The results of our calculations have been summarized 
in the Figure [2] The FT of the LDOS oscillations is 
plotted near a short-ranged diagonal impurity potential 
V(fc) ex L Given any other nontrivial form of this po- 
tential, the LDOS modulations may be found straight- 
forwardly from the above prescription. 

It is worth noting here that we have only quantified 
the oscillation parameters that may be observed in the 
spatial Fourier transform and not the Power Spectrum, 
examples of which are however also shown in Figure [2] 
(LI, Rl). Upon squaring the FT modulus to obtain the 
power spectrum the result could have twice as many os- 
cillations - this needs to be kept in mind when comparing 
the foregoing results with experimental signatures. 

In conclusion, in this work we have laid out the frame- 
work for calculating the LDOS modifications around an 
impurity in graphene in the presence of a strong mag- 
netic field. We use the linearly dispersing chiral quasi- 
particle theory. To calculate the QPI we have derived the 
graphene green's function in a magnetic field. There are 
two distinct regimes - in case of a strong field we have 
a situation of QHE while in the opposite case of a weak 



field the level broadening Y (due to lock-in modulation of 
STM voltage or due to impurity scattering, etc) will be 
larger then the Landau level splitting AE ~ B 1 ! 2 . We 
considered the case of a strong magnetic field. To this end 
we established a series expansion in Y / AE. In this limit 
our approach can be used to obtain the LDOS oscillations 
for any impurity potential. While the exact form of these 
oscillations vary by impurity type, we have identified a 
few important characteristics that may be observed in 
the FT of these oscillations — impurity-induced LDOS 
modulations in a magnetic field thus offers an alternative 
avenue for Landau level spectroscopy using local probes. 
(Note added: We recently became aware of the preprint 
[T2] where similar questions have been addressed.) 
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